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It usually writes the boundary condition of the wave equation in the Coulomb field as 
a rough form without considering the size of the atomic nucleus. The rough expression 
£ — \ brings on that the solutions of the Klein-Gordon equation and the Dirac equation with the 

. Coulomb potential are divergent at the origin of the coordinates, also the virtual energies, 

when the nuclear charges number Z > 137, meaning the original solutions do not satisfy the 
conditions for determining solution. Any divergences of the wave functions also imply that 
the probability density of the meson or the electron would rapidly increase when they are 
. closing to the atomic nucleus. What it predicts is not a truth that the atom in ground state 

would rapidly collapse to the neutron-like. We consider that the atomic nucleus has definite 
i' radius and write the exact boundary condition for the hydrogen and hydrogen-like atom, then 



00 
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newly solve the radial Dirac-Coulomb equation and obtain a new exact solution without any 
mathematical and physical difficulties. Unexpectedly, the K value constructed by Dirac 



33 . 

is naturally written in the barrier width or the equivalent radius of the atomic nucleus in 



solving the Dirac equation with the exact boundary condition, and it is independent of the 
quantum energy. Without any divergent wave function and the virtual energies, we obtain a 
new formula of the energy levels that is different from the Dirac formula of the energy levels 



in the Coulomb field. 
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I> . I. INTRODUCTION 
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The Dirac equation for the hydrogen atom has been treated in modern mathematical physics 
textbooks[l]. It is well known that the Dirac equation succeed in many respects. The Dirac 
equation is compatible with the theory of relativity, and it describes the spin of the electron and 
its magnetic moment in a completely natural way and so on. Especially, the distinguished Dirac 
formula of energy levels in Coulomb field can explain the fine-structure of the hydrogen atom. 
This is one of the important indicates of the achievements of the Driac theory. However, the Dirac 
operator is not apple-pie. According to the Ref.l, the main difficulty with a quantum mechanical 
on-particle interpretation is the occurrence of states with negative (kinetic) energy. Interaction 
may cause transitions to negative energy states, so that there is no hope for a stability of matter 
within that framework. On the other hand, we have to be aware of the fact that a quantum 
mechanical interpretation leads to inconsistencies if pushed too far. The localization problem and 
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the Klein paradox are still no clear solution, even in quantum electrodynamics. Investigating the 
Dirac equation, one should be not too far from Dirac's point of view: "... a book on the new 
physics, if not purely descriptive of experimental work, must be essentially mathematical". Here 
we start with the above reference and take new attention to one of the mathematical difficulties 
concealed in the Dirac equation with a Coulomb potential, which have not been obtained any 
logical mathematical explain, and we show finally that these paradoxes do not exists actually, such 
as the virtual energies and the divergence of all wave function on the ground state as r — ► in 
pure Coulomb field. We show that all of this kind of difficulties is due to the original incorrect 
mathematical methods for solving the radial Dirac equation. 

It is well known that, according to the boundary condition to solve Schrodinger equation [I]] 
for the hydrogen atom one can naturally obtain the Bohr formula[3] of energy levels, this is the 
prominent sign that the quantum mechanics is different from the classical mechanics. Considering 
the relativistic effect, Dirac introduced his relativistic wave equation for the single electron [4] 
in 1928. Darwin[5|] and GordonWj first obtain the exact solution of the Dirac equation with a 



Coulomb potential. Biedenharm 



3] , Wong and Yeh flfl, SuQ etc. also constructed the different 



ll(], Kalus and Wiist[12] 



second Dirac equation and obtain the different form of solutions. Nenciu 
investigated the different construction methods of self-adjoint extension of the Dirac operators with 
coulomb potential, and it is also showed that the distinguished self-adjoint extensions given by the 
two methods are identical. In history, it was ineffectual to use the Klein-Gordon equation [l^ \\\ ■ 
to describe the hydrogen atom because its eigenvalues of the quantum energy is incompletely agrees 
accurately with the experimentally observed hydrogen spectra. We investigated the mathematical 
foundation of constructing all wave equations for the quantum system by the numbers and found 
some mathematical problems that have been ignored in the Dirac equation. We know that all wave 
equations for the quantum system have not been strictly demonstrated in quantum mechanics, 
but the Schrodinger equation has not any mathematical difficulty, it is regarded as a fundamental 
assumption and is generally accepted. However the original solution of the Dirac equation in 
Coulomb field is divergent at the origin of coordinate and the energies for all atoms with nuclear 
charges number Z > 137 are the virtual numbers, all of these problems are actually serious mistakes 
from the incorrect mathematical methods. If a solution of some wave equation contains some 
mathematical or physical difficulties even mistakes, it should imply some new laws that have not 
been found \}\- The mathematical difficulties concealed in the Dirac equation with the Coulomb 
potential should imply some new conclusions that have not been found also. 

For briefness to discuss the essential of the problems, we only focus on the original exact solution 
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of the Dirac equation for the hydrogen atom in present paper. It is should be accentuated that the 

1 



boundary condition exerts decisive action in solving the wave equation^ 
However such important condition is not always attended in the history of the relativistic quantum 
mechanics That the whole Dirac wave functions of the hydrogen and 

hydrogen-like atom are divergent at the originof the coordinate implies that the atom on the 
ground state would collapse to the neutron-like. By all appearances, the deduction falls short of 
the true. In fact, the original mathematical methods for solving the radial Dirac equation are 
incorrect, although one can obtain the fine-structure formula of the energy-levels. This is why we 
have those mathematical difficulties. For finding the correct eigensolution of the some differential 
equations that cannot be optimized such as the radial Dirac equation with Coulomb, there is some 
theorems which need ones to prove 



31] 32]. Any wave function must accord with the conditions 
for determining solution and the physical signification. The virtual energies and the divergence of 
the original solution for the Dirac equation are due to the rough boundary condition. Writing the 
exact boundary condition and solving the Dirac equation outside the nucleus again, one can obtain 
a new exact wave function which has no any mathematical paradox, including the divergence of the 
wave function and the virtual energies. Howvere, the new strict mathematical deduction indicates 
that the K value constructed by Dirac is naturally written in the radius of the hydrogen atom, one 
can find that the new formula of the energy levels is no longer the Dirac formula, it is the correct 
result from correct mathematical method for solving the redial Dirac equation. 



II. ROUGH BOUNDARY CONDITION AND DIVERGENCE OF DRIAC FUNCTION 

Using the Schrodinger equation to describe the quantum system of the hydrogen atom, it usually 
does not consider the size of the atomic nucleus. The boundary condition of the atom with nuclear 
charges number Z is written as the rough form 

R(r ->0)^±oo, i?(r^oo)=0, -oo < R (0 < r < oo) < oo (1) 

where R is the radial wave function. Solving the Schrodinger equation by using this rough bound- 
ary condition, the Bohr formula of the energy levels is naturally obtained, it is one of the most 
consummate parts in quantum mechanics. In Dirac theory, because of absence for considering size 
of the atomic nucleus, one still uses the above rough boundary condition to solve the radial Dirac 
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equation and obtains the distinguished Dirac formula of the energy levels in the Coulomb field 33] 

„2 

^=^= (2) 



E 



mc 



1 + a 2 



n, 



+ VK 2 - Z 2 



a- 1 



where K = ±1, ±2, ±3, • • • , constructed by Driac, n r = 0, 1, 2, • • • , and a is the fine structure 
constant. It sticks out a mile, for the ground state or says S-state also, n = 0, K = ±1, when 
Z/,137, the energies of the system must be virtual numbers. This is pure mathematical problems. 
We don't think that those explain which departure mathematics too far are correct. 

On the other hand, the corresponding Dirac wave function takes the form with two components 
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where a = \lm 2 c^ — E 2 j he and the coefficients of the polynomial satisfy the system of the recursive 
relations 

\b v _ x + d„-x - Zab v -[K+ VK 2 - Z 2 a 2 + v^d v = 
b v -i + Xd v -i + (k~ ^K 2 - Z 2 a 2 -v)b u + Zad v = 
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where A 



mc l , % . However the wave function for the S state is divergent, as K = 1, whatever 



n r takes any value, the expression ([3]) becomes 
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this is a typical mathematical difficulty, which the solution balance out the condition for determin- 
ing solution of the wave equation. 

Any mathematical difficulty in physical theory must lead to some deductions that contravene 
the order of nature. When looking from a physical point of view, the divergence of the Dirac wave 
function for S state implies that the probability density of the electron around the nucleus rapidly 
increases as it close to the atomic nucleus. The probability density from the wave function with 
two components is defined as 



p (r, t) = R + (r, t) R (r, t) 



According to ([3]), we have 
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In this case the radial probability density of the electron for the relativistic hydrogen is as follows 
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For S-state which implies.fr = 1, as r — > 0, the above formula becomes 



lim p = oo 
r-»0 



(9) 



what this result predicts should be that the hydrogen and hydrogen-like atom in the ground state 
must rapidly collapse to the neutron-like. However the fact is not thusness. That is to say, the 
original solution of the Dirac equation for the hydrogen and hydrogen-like atom neither agrees 
with the mathematical principle nor agrees with the physical signification. Unexpectedly, such 



34] 35] 36] so that hardly might one open out 



divergence was defined as so-called "mild divergence" 
its actual meaning, and the correct deduction have been buried. We know that the Klein-Gordon for 
the meson without spin has the same divergence, but the Klein-Gordon divergence in the Coulomb 
field can be eliminated by the suitable mathematical method. Only one demonstrate some new 
theorems for finding the eigensolutions set of some differential equations with the variable coefficient 
can find the correct eigenvalues set of the corresponding wave equations (Ref.31 and Ref.32). 

Using a cut-off procedure for the potential that is similar to the case of considering an extended 
nucleus to blench the divergence should be independent of the exact solution for the Dirac equation 
with the Coulomb potential, and it oppresses the exact solution of the wave equation. For the exact 
solution, why coming forth the mathematical difficulty such as the expression ([5]) and ([9]) and the 
virtual energies is that the size of the nucleus of the hydrogen atom is not considered in the rough 
boundary condition (pQ), and the nuclear are regarded as the point in geometrical meaning. The 
point in geometrical meaning falls short of the actual case of the atomic nucleus. In fact, the 
necessary of the normalizable wave function of the hydrogen atom has been discussed home and 
widely in some modern physics textbookf^]. Now one should consider the actual size of the atomic 
nucleus to rewrite the boundary condition then find the eigensolution of the Dirac equation for the 
hydrogen and hydrogen-like atom. 



III. EXACT BOUNDARY CONDITION AND NEW SOLUTION OF DIRAC EQUATION 

Boundary conditions for any wave equations are written out basing on the structure of the 
physical model and distributing character of the physical quantity. Consider two basic facts, one 
is that the atomic nucleus has definite size, we suppose its equivalent radius or barrier width is 8. 
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Another is that the electron does not enter the inside of the atomic nucleus, and does not collide 
to and rub with the atomic nucleus. In this way, any wave equation that describes the atom has 
the same exact boundary condition 

R (r < S) ^ ±00, R (r -> 00) = 0, -00 < R {5 < r < 00) < 00 (10) 

one would recover the Bohr formula of the energy levels if uses this exact condition to solve the 
Schrodinger equation for the hydrogen atom. It is well known that the Schrodinger equation is 
very consummate in mathematics. However the new formula of the energy levels coming from the 
Dirac equation with the exact boundary condition is not as exact as the Dirac formula with the 
rough boundary condition and the divergence for the wave function. Of course, one would obtain 
the satisfying formula that is as exact as the distinguished Dirac formula when considering the 
spin-orbit coupling in the Dirac equation. It should be not one and only choice that for explaining 
the fine structure of the hydrogen atom one still uses the rough boundary condition for the Dirac 
equation and still blench those mathematical and physical difficulty such as the divergence and the 
virtual energies. One note that the boundary conditions are related to the self-adjointness of the 
operator, also one of methods to obviate the divergence of the Dirac function was given by Deck, 
Amar and Fralickj38f|. 

Now use the exact boundary condition (|10p to solve the radial Dirac equation of the hydrogen 
atom. It usually introduces a mathematical transformation 

/ F(r)/r \ 

R=[ U/ (11) 
\G(r)/r ) 

and translate the radial Dirac equation for the hydrogen atom 
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into the following form 
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where K = ±1, ±2, • • • Considering the exact boundary condition ()10p . introduce the transform 

£ = r - 5 (£ > 0) (14) 
the boundary (|10p can be overwritten as follows 

R (f 0) + ±00, R (£ — > 00) = 0, -00 < R (0 < £ < 00) < 00 (15) 
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then r = £ + 5, substituting it into (|T3[) . one obtains 



Tic + 5+5 ' r 
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Because .F— ►OjG— >0as£— ►oo, the system of differential equations (fTI)]) has the formal solutions 
with weight function of asymptotic solutions 



(17) 



Substituting (fTTj) into equation (|T6j) . one then obtains 
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/ + (£ + S)f + [#-a(£ + £)]<7 = 
9-(C + d)^ + [K + (a^ + 5)]f = 

the Eigensolutions of equations (|18p correspond to quantum energy are two interrupted serieswhich 
the number of terms is determined by the eigenvalues. 

In order to find the general series solutions for equations (|18p . it is assumed that the formal 
solutions are 
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Substituting into equations (fTHj) . one obtains the linear system of recursive relations 
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hence the coefficient of the power series satisfy the following system of recursive relations 



(20) 



+5 (a + v + 1) c^+i + (If + a + - 5a) d v = 
^^_ 1 + (^<5 + a) d. + a^a 



(21) 



-5 (a + v + 1) 6^+1 + (if - a - v + 5a) b v = 

Corresponding to v = — 1 the indicial equations are given that oo"6o = and 5ado = 0. Because 
5 / 0, &o 7^ and do ^ 0, one obtains 



a = 



(22) 
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so that the wave functions satisfy the boundary condition at r — ► 5 namely £ — ► 0, the above 
equations reduce to 

M =^-K-\ + (^j^<J + a) b v - ad u -! + <5 (u + 1) + (if + i/ - oa) d v = 
+ (^^5 + ajd u + ab,,-! -S(y + 1) b u+1 + {K-u + 5a)K = 

Respectively evaluate for v = 0, 1, 2, • ■ ■ , n r , b nr+ \ = d nr+ \ = 0, make use of that 6_2 = d_2 = 
and b-x = d_i = 0, equations (|23|) give 

^5 + a) 6 + (K - 5a) d + Sdx = 
(if + 5a) & + f M ¥^6 + a) do - = 



^ he 

M =^b Q + (^f^<5 + ajh-ad + (K + l- 5a) d x + 25d 2 = 
ab + (K-l + 5a) h + ^^d + ( M± ^ L 8 + a) d x - 25b 2 = 



(24) 



E—mc 2 i i | E—mc 2 



Hc— b nr-l + + a)b nr - ad nr -i + (K + n r - 5a) d Ur = 



a6 nr _! + {K-n r + 5a) b nr + ^p^-i + (^£^5 + a) d Ur = 
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The last two formulas are linearly dependent. Use (E + m.c 2 ) /Tic to multiply the third formula 
from bottom and use \/m. 2 c 4 — E 2 I he to multiply the fourth formula from bottom, and then add 
the two new formulas, it is given as follows 



+ 



a (E + mc 2 ) + (K - n r ) Vm 2 c 4 - E 2 
(K + n r ) (E + mc 2 ) + a\/m 2 c 4 - E 2 



(25) 

d n „ = 



Substituting for the second formula from the bottom, one will obtain a new formula of the energy 
levels for the hydrogen atom 

2 

TYIC 

E= ; = (n r = 1,2,3,---) (26) 



^ 1 ' ''I' 

One can strictly demonstrate that n r > 1. It is different from the Dirac formula of the energy 
levels for the hydrogen atom. This result is the inevitable deduction of the Dirac equation with the 
exact boundary condition for the hydrogen atom. With the exact boundary condition (|10p and the 
new formula of the energy levels (|26p . all of the corresponding wave functions satisfy the boundary 
conditions and there is not any virtual energies. 
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According to (jTTj) . (jHj), (|T6j) . (|19j) . (j23|) . the whole wave function with the exact boundary 
condition is as follows 



R 



( e a i t t!/ \ 



i/=i 



5 £ <W* 



, (n r > 1) (27) 



the coefficients of corresponding polynomial are determined by the system of recursive relations 
(|23p . All appearance, at the boundary of the hydrogen atom 

/ Constant \ ( 0\ 

lim R = , lim R = (28) 

^° y Constant J *-"*> \ J 

Make use of the definition ([6]), the probability density of the electron appearing outside the nucleus 
of the hydrogen atom takes the form 

/ p -aZ nr \ 2 ( p -a£ n r \ 2 



i/=0 / \ i/=0 



homoplastically, one has 



lim p = Constant, lim p = (30) 
IV. CONCLUSIONS 

In this paper we expatiated that the divergences of the Dirac function and the virtual energies 
of the Dirac formula of the energy levels for the hydrogen and hydrogen-like atom are due to the 
traditional rough boundary condition. By using the exact boundary condition one can obtain a new 
solution of the Dirac equation in the Coulomb field. The new solution without any mathematical 
difficulty gives a new formula of the energy levels which is different from the distinguished Dirac 
formula. One can find that in the new solution the K values constructed by Dirac have been 
written in the radius of the atomic nucleus and they are independent of the new formula of the 
energy levels. Only when looking from point of view for the formula of the energy levels, the new 
formula is not as exact as the Driac formula. However, considering the spin-orbit coupling or some 
new potential parameters, one will obtain the exacter formula. 

It is well known that the Dirac equation succeed in many important prognostics 41]. Dirac's 



equation has an infinite number of solutions with negative energies. It leads to discover the existence 
of the positron and implies an unexpected relativistic interaction between an electron's translational 
motion and spin, which leads to a violent oscillation of the particle at very high frequencies and 
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over distances of roughly one Compton wavelength. For Dirac's electrons, zitterbewegung takes 
place whenever the electron wave function includes both positive and negative energy components. 
It takes both sets of states to build up an arbitrary electronic state. The distillate of quantum 
mechanics is to naturally obtain the formula of the energy levels for the bound state by solving the 
wave equation with the boundary condition. Disclosing the mathematical and physical difficulty 
concealing in the Dirac equation with the Coulomb potential does not imply to negative the Dirac 
theory. Our works wishes to vindicate the mathematical logic for the wave equation very precise, 
at least to eliminate the divergence by using the exact boundary condition to replace the rough 
boundary condition. 

No matter which wave equation is used to describe the quantum system of the bound statge, their 
eigenvalues set and eigensolutions set must be in agreement with the uniqueness, and the solution 
must accord with the conditions for the exact solution but not any approximate solution such as cut 
off potential. In principle, we cannot immolate the mathematical rule to obtain a formula only for 
agreeing with the experimentally observed hydrogen spectra. Actually the divergence of the original 
solution for the Dirac equation in the Coulomb field is nonexistent. The Dirac equation is more 
and more widely applied for various models. One should note that using different mathematical 
methods and different boundary condition to solve the differential equation will obtain different 
results 42 



431 ] . The rough boundary condition for the Dirac equation with the Coulomb potential 
brings on so mathematical contradictions. We need to revise all incorrect mathematical methods 
and search the correct mathematical methods to find the correct and exact solution of the various 
wave equations. One can find that the classical solution of the plane transverse electromagnetic 
mode of the Maxwell equation is also incorrect 44j] , because of the incorrect mathematical method 



for solving the wave equation. For quantum mechanics, only the exact solution of the wave equation 
can we know if the differential wave equation has the eigenvalues connecting the energy levels. 
Further research should explore the exacter formula of the energy levels for the hydrogen and 
hydrogen-like atom by using the exact boundary condition that accords with the structure of the 
atoms. Here it should be pointed out that the Schrodinger equation and the Klein-Gordon equation 
in the Coulomb field with the exact boundary condition have the same corresponding formula of 
the energy levels. 
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